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Abstract 

Let F := (/i, . . . , f n ) be a random polynomial system with fixed n-tuple of supports. 
Our main result is an upper bound on the probability that the condition number of / 
in a region U is larger than 1/e. The bound depends on an integral of a differential 
form on a toric manifold and admits a simple explicit upper bound when the Newton 
polytopes (and underlying covariances) are all identical. 

We also consider polynomials with real coefficients and give bounds for the expected 
number of real roots and (restricted) condition number. Using a Kahler geometric 
framework throughout, we also express the expected number of roots of / inside a 
region U as the integral over U of a certain mixed volume form, thus recovering the 
classical mixed volume when U = (C*) n . 
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1 Introduction 

From the point of view of numerical analysis, it is not only the number of complex solutions 
of a polynomial system which make it hard to solve numerically but the sensitivity of its 
roots to small perturbations in the coefficients. This is formalized in the condition number, 
M/j (°f- Definition [| of Section |1 . 1| ) , which dates back to work of Alan Turing ||Tur36| . 



In essence, fi(f, Q measures the sensitivity of a solution ( to perturbations in a problem 
/, and a large condition number is meant to imply that / is intrinsically hard to solve 
numerically. Such analysis of numerical conditioning, while having been applied for decades 
in numerical linear algebra (see, e.g., [ Dem97|| ), has only been applied to computational 



algebraic geometry toward the end of the twentieth century (see, e.g., [|SS93b|| ) . 

Here we use Kahler geometry to analyze the numerical conditioning of sparse polynomial 
systems, thus setting the stage for more realistic complexity bounds for the numerical solu- 



tion of polynomial systems. Our bounds generalize some earlier results of Kostlan |Kos93 
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and Shub and Smale ||SS96|| on the more restricted dense case, and also yield new formulae for 



the expected number of roots (real and complex) in a region. The appellations "sparse" and 
"dense" respectively refer to either (a) taking into account the underlying monomial term 
structure or (b) ignoring this finer structure and simply working with degrees of polynomi- 
als. Since many polynomial systems occuring in practice have rather restricted monomial 
term structure, sparsity is an important consideration and we therefore strive to state our 
complexity bounds in terms of this refined information. 

To give the flavor of our results, let us first make some necessary definitions. We must first 
formalize the spaces of polynomial systems we work with and how we measure perturbations 
in the spaces of problems and solutions. 

Definition 1 Given any finite subset Ac Z n , let ^(A) (resp. J-"-r(A)) denote the vector 
space of all polynomials in C[x%, . . . , x n ] (resp. WL[x\, . . . , x n }) of the form Yl c a x a where the 

notation x a := x ai ■ ■ ■ x an is understood. For any finite subsets A\, . . . , A n C Z n we then let 
A := (Ax,..., A n ) and F C {A) := x • • • x ^x:(A») ( re $P- ^k(-A) := ^k(Ai) x • • • x 

F R (A n )). o 

The n-tuple A will thus govern our notion of sparsity as well as the perturbations allowed 
in the coefficients of our polynomial systems. It is then easy to speak of random polynomial 
systems and the distance to the nearest degenerate system. Recall that a degenerate root 
of / is simply a root of / having Jacobian of rank < n. 

Definition 2 By a complex (resp. real) random sparse polynomial system we will 
mean a choice of A := (A\, . . . , A n ) and an assignment of a probability measure to each JFc(Aj) 
(resp. J~m.(AA) as follows: endow JF c (y4j) (resp. F^Ai)) with an independent complex (resp. 
real) Gaussian distribution having mean O and a (positive definite and diagonal) covariance 
matrix Ci. Finally, let the discriminant variety, T,(A), denote the set of all f G Tc(A) 
(resp. fEJ r ^(A)) with a degenerate root and define J-'^A) := {f & T<c(A) \ f(() — 0} (resp. 
F ( {A):={feMA)\f{Q = 0})andZ ( {A):=F ( {A)nZ{A). o 

Theorem 1 Suppose Acl/ 1 is a finite set with a convex hull of positive volume and A: = 
(A, . . . ,A). Then there is a natural metric d(-, •) on J-"c(A) such that //(/, Q = ^j^^yy- 

n 

Furthermore, 



Prob 



KfX)>\for some root (e(C*) n off 



< n 3 (n + l)Vol(A)(#A - l)(#A - 2)e 4 



where f is a complex random sparse polynomial system, j^A denotes the number of points in 
A, and Vol(A) denotes the volume of the convex hull of A (normalized so that Vol(0, e\, . . . , e n ) ■ 

V- 

The above theorem is in fact a simple corollary of two much more general theorems (Theorems 
H and [5D which also include as a special case an analogous result of Shub and Smale in the 



dense case [|BCSS98| , Thm. 1, Pg. 237]. We also note that theorems such as the one above 



are natural precursors to explicit bounds on the number of steps required for a homotopy 



2 



algorithm 



to solve /. We will pursue the latter topic in a future paper. Indeed, 



one of our long term goals is to provide a rigourous and explicit complexity analysis of the 
numerical homotopy algorithms for sparse polynomial systems developed by Verschelde et. 
al. |WC94| , Huber and Sturmfels flHS9lj and Li and Li jX0l| . 

The framework underlying our first main theorem involves Kahler geometry, which is 
the intersection of Riemannian metrics and symplectic and complex structures on manifolds. 
On a more concrete level, we can give new formulae for the expected number of roots of / 
in a region U. For technical reasons, we will mainly work with logarithmic coordinates. 
That is, we will let T n be the n-fold product of cylinders (R x (R mod 2-7r)) n C C n , and 
use coordinates p + iq :— (pi + iqi, . . . , p n + iq n ) G T" to stand for a root ( : = exp(jo + iq) : = 
(e pi+m , . . . , e Pn+tqn ) of /. Roots with zero coordinates can be handled by then working in a 
suitable toric compactification and this is made precise in Section 0. The idea of working 
with roots of polynomial systems in logarithmic coordinates seems to be extremely classical, 
yet it gives rise to interesting and surprising connections (see the discussions in | |MZa| , MZb, 
[VrrOOH ). 

Theorem 2 Let A\, . . . , A n be finite subsets of7L n and UC^T n be a measurable region. Pick 
positive definite diagonal covariance matrices C\, . . . ,C n and consider a complex random 
polynomial system specified by the data (Al, Ci, . . . , A n , C„). Then there are natural real 
2-forms lu Ai , ■ ■ ■ , u An 071 T n such that the expected number of roots of f in exp U C (C*) n is 
exactly 

/_-j\n(n-l)/2 



7T 



UJ Al A • • • A U An 



In particular, when U = (C*) n , the above expression is exactly the mixed volume of the convex 
hulls of Ai, . . . , A n (normalized so that the mixed volume of n standard n-simplices is 1). 

See | BZ88| , SY93 ] for the classical definition of mixed volume and its main properties. The 
result above generalizes the famous connection between root counting and mixed volumes 



discovered by David N. Bernshtein [Ber76|. The special case of unmixed systems with iden- 



tical coefficient distributions (A\ 



A n , C\ 



C n ) recovers a particular case 



of Theorem 8.1 in |[EK95|| . However, comparing Theorem |2| and |[EK95| , Theorem 8.1], this 
is the only overlap since neither theorem generalizes the other. The very last assertion of 
Theorem |2] (for uniform variance Ci = I for all i) was certainly known to Gromov |Gro' 



and a version of Theorem ^| was known to Kazarnovskii [|Kaz81| , p. 351] and Khovanskii 
[Kho91| , Prop. 1, Sec. 1.13]. In | |Kaz81| |, the supports Ai are even allowed to have complex 



exponents. However, uniform variance is again assumed. His method may imply this special 
case of Theorem^, but the indications given in [ Kaz81 | were insufficient for us to reconstruct 
a proof. 

As a consequence of our last result, we can also give a coarse estimate on the expected 
number of real roots in a region. 

Theorem 3 Let U be a measurable subset ofM. n with Lebesgue volume \{U). Then, following 
the notation above, suppose instead that f is a real random polynomial system. Then the 
average number of real roots of f in exp U C WL is bounded above by 



(4 7 r 2 )- n/2 v / A(f7y- 



' (p,q)eUx[0,2Tv) r ' 



(-l)n(n-W UAl A ••• AiU An . 
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This bound is of interest when n and U are fixed, in which case the expected number of 
positive real roots grows as the square root of the mixed volume. 



1.1 Stronger Results Via Mixed Metrics 

Our remaining new results, which further sharpen the preceding bounds and formulae, will 
require some additional notation. 

Definition 3 We define a norm on jF c (y4j) by ||/ l || c -i := f l C~ l (f l ) H where, in the last 

expression, we consider f % via its row vector of coefficients and (-) H denotes the usual Her- 
ri 

mitian conjugate transpose. Finally, we define a norm on J-"c(A) by \\f\\ := WFWr- 1 ! an d 

1 i 
1=1 

a metric dn> on the product of projective spaces ¥(J-'c{A)) : = P(jF c (yli)) x • • • x ¥(J r c(A n )) by 

n II/*"- Xg^ 1 1 

dw(f, g) '-=Y1 m * n \\p\\ > where we implicitly use the natural embedding of ¥(T<c(Ai)) into 
i=i a gc* 11 

the unit hemisphere of Tc{Aj). o 

Each of the terms in the sum above corresponds to the square of the sine of the Fubini 
(or angular) distance between f l and g l . Therefore, dp is never larger than the Hermitian 
distance between points in T<c{A), but is a correct first-order approximation of the distance 
when g -> / in P(J r c (^)) (compare with [ |BCSS98| , Ch. 12]). 



Recall that T P M denotes the tangent space at p of a manifold M. 

Definition 4 Define the evaluation map, ev A , as follows: 

ev A : T x T n -> C n 

((/ 1 ,--->/ n )>P + *9) (/ 1 ( ex P(P + ^)) ) --- ; / n ( ex P(P + ^)))- 

Given any root exp(p + iq) of an f in J-"c(A), the condition number of / at p + iq, 
fx(f,p + iq), is then defined to be the operator norm 



DG\ 



max 

yi=i 



DG 



where G is the unique branch of the implicit function which satisfies G(f) = p + iq and 
evji(g,G(g)) = O for all g sufficiently near f, and DG : T/jF c (yi) — > T p+iq T n is the 
derivative of G. (We set the condition number fi(f,p + iq) := +oo in the event that Df is 
zero and G thus fails to be uniquely defined.) o 

Note that the implied norm on TfT<c{A) was detailed in the previous definition, while the 
implied norm on T p+ i q T n has intentionally been left unspecified. This is because while 
J-"c(A) admits a natural Hermitian structure, the solution-space T n admits n different nat- 
ural Hermitian structures (one from each support Aj, as we shall see in the next section). 
Nevertheless, we can give useful bounds on the condition number and give an unamibiguous 
definition in certain cases. 

Theorem 4 (Condition Number Theorem) If(p,q)£ T n is a non- degenerate root of f 
then 

max min \\DG f f\\ Al < , „ r < max max \\DG f f\\ Ai ■ 

II/II<1 1 «PU>^(p,g)J ||/||<l 1 



In particular, if A\ = ■ ■ ■ = A n and G\ = ■ ■ ■ = C n , then 

1 



max min H^DG//!!^ = max max H-DC//!^ 



[|<1 ' ' * |[/||<1 ' Mf^( P ,q)) 

and we can define fi(f', (p, q)) to be any of the three preceding quantities. 

This generalizes PCSS98| , Thm. 3, pg. 234] which is essentially equivalent to the last assertion 
above, in the special case where is an n-column matrix whose rows {^4"} a consist of all 



partitions of cL into n non-negative integers and Cj = Diag„ ■ - , , „,, — , ,,, , — ^ — - . ... ,, ■ 

- in short, the case where one considers complex random polynomial systems with f % a 
degree di polynomial and the underlying probability measure is invariant under a natural 
action of the unitary group U(n+1) on the space of roots. The last assertion of Theorem ^ 
also bears some similarity to Theorem D of [ Ped96|| where the notion of metric is considerably 
loosened to give a statement which applies to an even more general class of equations. 
However, our philosophy is radically different: we consider the inner product in Tc(A) as 
the starting point of our investigation and we do not change the metric in the fiber T( p , q )- 
Theorem |4] thus gives us some insight about reasonable intrinsic metric structures on T". 

In view of the preceding theorem, we can define a restricted condition number with 
respect to any measurable sub-region U C T n as follows: 

Definition 5 We let u(f; U) := — \ r. Also, via the natural GL(n)-action on 

T^ q )T n defined by (p, q) \— > (Lp, Lq) for any L e GL{n), we define the mixed dilation of the 
tuple (u Al , ■ ■ ■ ,uJ An ) as: 



k(w Ai , ■ ■ ■ , u An ; (p, q)) := min max 



max\\ u \\=i(u) Ai )(p jq )(Lu, JLu) 



L€GL(n) i min M=1 (uj Ai ) M (Lu, JLu) 

where J : TT n — > TT n is canonical complex structure of T n . Finally, we define Ku : = 
sup(p 9 ) eC/ k(uj Ai , ■ ■ ■ ,u An ; (p,q)), provided the supremum exists, and K[/: = +oo otherwise. 

We can then bound the expected number of roots with condition number p > on 
U in terms of the mixed volume form, the mixed dilation Kjj and the expected number of 
ill-conditioned roots in the linear case. The linear case corresponds to the point sets and 
covariance matrices below: 



.4 



Lin 





1 



a 



Lin 



Theorem 5 (Expected Value of the Condition Number) Let v (n, e) be the proba- 
bility that a complex random system of n polynomial in n variables has condition number 
larger than e^ 1 . Let v A {U,e) be the probability that fj,(f,U) > e^ 1 for a complex random 



polynomial system f with supports A±, 



Then, v A (U,e) < J^" A 



A n and covariance matrices C\, 



^ A Lin 



Kye). 
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Our final main result concerns the distribution of the real roots of a real random polyno- 
mial system. Let i^(n, e) be the probability that a real random linear system of n polynomials 
in n variables has condition number larger than e . 

Theorem 6 Let A — A\ — ■ ■ ■ — A n and C = C\ — ■ ■ ■ = C n and let U CR™ be measurable. 
Let f be a real random polynomial system. Then, 

Prob [fi(f,U) > e- 1 ] < E(U) i*(n,e) 

where E(U) is the expected number of real roots on U . 

Note that E(U) depends on C, so even if we make U = M n we may still obtain a bound 
depending on C. Shub and Smale showed in | fSS93a| | that the expected number of real roots 
in the dense case (with a particular choice of probability measure) is exactly the square root 
of the expected number of complex roots. The sparse analogue of this result seems hard 
prove even in the general unmixed case: Explicit formula? for the unmixed case are known 



only in certain special cases, e.g., certain systems of bounded multi-degree [[Roj96| , [McLOO 



Hence our last theorem can be interpreted as another step toward a fuller generalization. 



2 Symplectic Geometry and Polynomial Systems 
2.1 Some Basic Definitions and Examples 

For the standard definitions and properties of symplectic structures, complex structures, 
Riemannian manifolds, and Kahler manifolds, we refer the reader to ||MS98| , |CCL99|| . A 



treatment focusing on toric manifolds can be found in ||Gui94| , Appendix A]. We briefly 



review a few of the basics before moving on to the proofs of our theorems. 

Definition 6 (Kahler manifolds) Let M be a complex manifold, with complex structure 
J and a strictly positive symplectic (1, 1) -form u on M (considered as a real manifold). We 
then call the triple (M,u, J) a Kahler manifold, o 

Example 1 (AfRne Space) We identify C M with M? M and use coordinates Z l = X' 1 + 
y/—lY % . The canonical 2-form ujz = X]f=i A dYi makes C M into a symplectic manifold. 

The natural complex structure J is just the multiplication by \f—\. The triple (C M , Uz, J) 
is a Kahler manifold, o 

Example 2 (Projective Space) Projective space P A/_1 admits a canonical 2-form defined 
as follows. Let Z = (Z\ ■ ■ ■ , Z M ) e (C M )*, and let [Z] = (Z 1 : • • • : Z M ) e P^ 1 be the 
corresponding point in P A/_1 . The tangent space T[^]P A/_1 may be modelled by Z 1 - C T^C M . 
Then we can define a two-form on p^^ 1 by setting: 

LO[ Z ](u,v) = \\Z\\~ 2 ujz(u, v) , 

where it is assumed that u and v are orthogonal to Z . The latter assumption tends to be 
quite inconvenient, and most people prefer to pull U[z] back to C M by the canonical projection 
7i : Z i— > [Z] . It is standard to write the pull-back r = tt*u)\z\ as: 

T z = ~^dJ*d ilog||Z|| 2 , 
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using the notation dr] = ^ ^ A dpi + ^ A dqi, and where J* denotes the pull-back by J. 

Projective space also inherits the complex structure from C M . Then U[z] is a strictly 
positive (1, l)-form. The corresponding metric is called Fubini-Study metric in C M or <£ M ~ l . 

o 

Remark 1 Some authors prefer to write y/ —ldd instead of — \dJ*d. The following notation 
: i %AdZ l andd V = Z l % 

7 / £. dZi A dZ, ZidZi A J2j ZjdZ, 



is assumed: drj = ~z. A dZi and drj = ^ ^ A dZi. Then they write tz as: 



Example 3 (Toric Kahler Manifolds from Point Sets) Let A be any M x n matrix 
with integer entries whose row vectors have n-dimensional convex hull and let C be any 
diagonal positive definite n timesn matrix. Define the map Va from C n into C M by Va '■ 



C l/2 



A 1 "1 



AM 

z A 



We can also compose with the projection into projective space to obtain 



a slightly different map Va = ft ° Va '■ C n — > P M_1 defined by Va '■ z i— > When C 

is the identity, the Zariski closure of the image of Va is called the Veronese variety and the 
map Va is called the Veronese embedding. Note that Va is not defined for certain values of 
z, like z — 0. Those values comprise the exceptional set which is a subset of the coordinate 
hyperplanes. 

There is then a natural symplectic structure on the closure of the image of Va, given 
by the restriction of the Fubini-Study 2- form r: We will see below (Lemma \3\) that by our 
assumption on the convex hull of the rows of A, we have that DVa is of rank n for z e (C*)" . 
Thus, by the above lemma, we can pull-back this structure to (C*) n by Qa — VX T - Also, we 
can pull back the complex structure o/P M_1 ; so that Qa becomes a strictly positive (1,1)- 
form. Therefore, the matrix A defines a Kahler manifold ((C*) K ,Qa, J)- o 

Lemma 1 Let A be a matrix with non-negative integer entries, such that Coirv(A) has 
dimension n. Then (Dva) p is injective, for all p e M n . 

Proof: The conclusion of this Lemma can fail only if there are pel" and ii / with 
(Dva) p u = 0. This means that 

P v A ( p) diag(w a) P Au = . 

This can only happen if di&g(vA) p Au is in the space spanned by (va) p , or, equivalently, 
Au is in the space spanned by (1, 1, • • • , 1) T . This means that all the rows a of A satisfy 
au = A for some A. Interpreting a row of A as a vertex of Conv(A), this means that Conv(A) 
is contained in the affine plane {a : au = A}. ■ 

Example 4 (Toric Manifolds in Logarithmic Coordinates) For any matrix A as in 
the previous example, we can pull-back the Kahler structure of ((C*) n , Qa, J) to obtain an- 
other Kahler manifold (T n , uja-, J)- (Actually, it is the same object in logarithmic coordinates, 
minus points at "infinity" '.) An equivalent definition is to pull back the Kahler structure of 

the Veronese variety by %a = Va ° exp. o 
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Remark 2 The Fubini- Study metric on C M was constructed by applying the operator — \d J* d 
to a certain convex function (in our case, | log H^H 2 ^ This is a general standard way to con- 
struct Kahler structures. In ^Gro9Qj , it is explained how to associate a (non-unique) convex 
function to any convex body, thus producing an associated Kahler metric, o 

For the record, we state explicit formulae for several of the invariants associated to the 
Kahler manifold (T n ,u>A, J)- First of all, the function g^ = g ° va is precisely: 



Formula 2.1.1: The canonical Integral gA (or Kahler potential) of the convex set asso- 
ciated to A 

9a{p) ■= g lo S (( ex P(^- • P)) T c (exp(A • p))J 

The terminology integral is borrowed from mechanics, and it refers to the invariance of 
gA under a [0, 2II) n -action. Also, the gradient of gA is called the momentum map. Recall 
that the Veronese embedding takes values in projective space. We will use the following 
notation: Va{p) = va(p) I\\va{p)\\- This is independent of the representative of equivalence 
class va(j>). Now, let va(p) 2 mean coordinatewise squaring and va(p) 2T be the transpose of 
va(p) 2 - The gradient of gA is then: 



Formula 2.1.2: The Momentum Map associated to A 

Vg A = v A (p) 2T A 



Formula 2.1.3: Second derivative of g& 

D 2 g A = 2Dv A {p) T Dv A {p) 
We also have the following formulae: 



Formula 2.1.4: The symplectic 2-form associated to A: 

(wa)(p, 9 ) = g ^2( D2 9A)ijdpi A dqj 



Formula 2.1.5: Hermitian structure of T" associated to A: 

((u,w) A ) M = u H (^D 2 g A ) P w 

In general, the function va goes from T n into projective space. Therefore, its derivative 
is a mapping 
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{Dv Ai 



. T t-r'n ( rji 
. l( n „\l — > 1, 



ip,q) • 1 (p,q) 



P 



M-l 



va(p + q^ 



■IV c 



For convenience, we will write this derivative as a mapping into C M , with range va{p + 
q\^l) ± . Let P v be the projection operator 



P, 



tVV 



II 



We then have the following formula. 



Formula 2.1.6: Derivative of va 



{Dva, 



(p,q) 1 VA(p+qv-l) 



T)Diag 



\vA(p + q^ 



A 



An immediate consequence of Formula 2.1.6] is: 



Lemma 2 Let f G J- a and (p,q) G T n be such that f ■ Va{p + gv^~T) = 0. T/ien, / ■ 



In other words, when (f oexp)(p+q\ / — 1) vanishes, -Df a and -Dva are the same up to scal- 
ing. Noting that the Hermitian metric can be written ((u, w) A )(p iq ) = u h DvA{p, q) H Dva(p, q)w, 
we also obtain the following formula. 



Formula 2.1.7: Volume element of (T n ,u>A, J) 

dT2 = det ( - D 2 g A {p) j dp x A • • • A dp n A dq x A ■ ■ ■ A dq n 



2.2 Toric Actions and the Momentum Map 

The momentum map, also called moment map, was introduced in its modern formulation by 
Smale [[5ma70|l and Souriau [|Sou70|l . The reader may consult one of the many textbooks in 
the subject (such as Abraham and Marsden [ |AM78|| or McDuff and Salamon ||MS98|| ) for a 
general exposition. 



In this section we instead follow the point of view of Gromov | Gro90 |. The main results 
in this section are the two propositions below. 

Proposition 1 The momentum map VgA maps T n onto the interior of Conv (A). When 
VgA is restricted to the real n-plane [q = 0] C T n , this mapping is a bisection. 

This would appear to be a particular case of the Atiyah-Guillemin-Sternberg theorem 
Ati82| , |GS82|| . However, technical difficulties prevent us from directly applying this result 



here.F 



lr The Atiyah-Guillcmin-Stcrnberg applies to compact symplectic manifolds and the implied compactifica- 
tion of T n may have singularities. 
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Proposition 2 The momentum map is a volume-preserving map from the manifold 

(T n ,ujA, J) into Conv(A), up to a constant, in the following sense: if U is a measurable 
region of Cony (A), then 

Vol((Vg A )- 1 (U))=ir n Vol(U) . 
Proof of Proposition [21: Consider the mapping 

|Conv(A) x T n 



M : T" 

(p,q) 



!|V0a(p),?) 



Since we assume dimConv(A) = n, we can apply Proposition |T] and conclude that M is 
a diffeomorphism. 

The pull-back of the canonical symplectic structure in R 2n by M is precisely u A , because 
of Formulae [2.1.3| and |2.1.4j . Diffeomorphisms with that property are called symplectomor- 
phisms. Since the volume form of a symplectic manifold depends only of the canonical 
2-form, symplectomorphisms preserve volume. We compose with a scaling by \ in the first 
n variables, that divides Vol(£7) by 2 ra , and we are done. ■ 

Before proving Proposition [I], we will need the following result about convexity which has 
been attributed to Legendre. (See also ||Cro90| , Convexity Theorem 1.2] and a generalization 
in ||Avi76|, Th. 5.1].) 



Legendre's Theorem If f is convex and of class C 2 on M. n , then the closure of the image 
{Vf r : r G 1"} in M™ is convex. 

By replacing / by g A , we conclude that the image of the momentum map VgA is convex. 
Proof of Proposition |]: The momentum map V#a maps T n onto the interior of Conv(A). 
Indeed, let a = A a be a row of A, associated to a vertex of Conv(A). Then there is a direction 



v e 



such that 



a - v = max x ■ v 

xeConv(A) 



for some unique a. 

We claim that a e V5u(B^ n )- Indeed, let x(t) = VA{tv), t a real parameter. If b is another 
row of A, 



e ta-v > e tb-v 



as t — > oo. We can then write VA{tv) 2T as: 







T 




v A {tv) = 




CDiag 





Since C is positive definite, C aa > and 

v A {tv) 2T 



lim VA(tv 



,2T 



lim , 



T ^ act T 
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where e a is the unit vector in IR M corresponding to the row a. It follows that lim^oo V(yu(tt>) = 
a 

When we set q = 0, we have det D 2 gA ^ on M. n , so we have a local diffeomorphism at 
each point p G M n . Assume that (V(?a) p = (V(yu)p' for p ^ p'. Then, let j(t) = (1 —t)p + tp'. 
The function t i— ► (V^) 7 m7'(t) has the same value at and at 1, hence by Rolle's Theorem 
its derivative must vanish at some t* G (0, 1). 

In that case, 

(^)#)(7'(n,7'(n) = o 

and since 7'(i*) = p' - p ^ 0, det-D 2 ^ must vanish in some p G M n . This contradicts 
Lemma [L| ■ 



2.3 The Condition Matrix 



Following |PCSS98[| , we look at the linearization of the implicit function p + qy/—l = G(f) 
for the equation evji(f,P + qV~ T) = 0. 

Definition 7 T/te condition matrix of ev at (f,p + q^—l) is 

DG = D T n(ev)~ l D T (ev) , 

where T = Ta x x • • • x jF An . 

Above, D T n[ev) is a linear operator from an n-dimensional complex space into C n , while 
D T (ev) goes from an (Mi + • ■ ■ + M n )-dimensional complex space into C n . 

Lemma 3 If p + iqE T n and /(exp(p + zg)) = O t/ien 

det (DGDG H y l d Pl Ad qi A---A dp n A dq n = (-l)^™- 1 )/ 2 /\ y/^if . (Dv Ai ) M dpA 

A f ■ {Dv At ) {v ,- q) dq . 

Note that although /' • (DvA^^dp is a complex-valued form, each wedge f l ■ (DvAi)( P , q )dp A 
f l ■ {DvA i )[ P - q )dq is a real-valued 2-form. 
Proof of Lemma |3|: We compute: 



DJev 



\(p,q) 



and hence 
Also, 



■££i^(p + ?V=T)4E 

E*:A(p + q^l)dfS\ 
D r {ev)D r {ev) H = diag ||£aJ| 2 - 



D T n(ev) 



7 1 • Dva 1 

r ■ Dv An 
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Therefore, 



det 



det (DG^DGlJ- 1 
We can now use Lemma § to conclude the following: 



ll'vijl ' 



Formula 2.3.1: Determinant of the Condition Matrix 



det (DGmDG*J 



det 



f l -Dv Al 
f n ■ Dv An 



We can now write the same formula as a determinant of a block matrix: 

7 1 • Dv Al 



det (DG M DG» q) ) * = det 



f n ■ Dv Ar> 



f 1 ■ Dv Al 



f n ■ Dv An , 



and replace the determinant by a wedge. The factor (— 1)"(" -1 )/ 2 comes from replacing 
dpi A • • • A dp n A dqi A • • • A dq n by dpi A dqi A • • • A dp n A dq n . ■ 
We are now ready to prove our main theorems. 



3 The Proofs of Theorems 1—6 

We will prove our main theorems in the following order: 1, 2, 4, 5, 3, 6. 
3.1 The Proof of Theorem 0] 

The first assertion, modulo an exponential change of coordinates and using the multi- 
pro jective metric dp(-, •), follows immediately from Theorem^. 

As for the rest of Theorem [1], Theorem |] applied to the linear case then provides the 
following interpretation of u Lin (n, e): 

z/ Li >,e) = Prob [d P (/,E (p ,, ) )<e], 

where / is a complex random linear polynomial system, and (p, q) is such that /(exp(p + 
iq)) = 0. So we are on our way to proving the inequality 

Prob [d ¥ (f, E M ) < e] <n\n + l)Vol(A)(#A - - 2)e 4 , 
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for general /, which clearly implies our desired bound. 

To prove the latter inequality, recall that by the definition of the multi- projective distance 
rfp(-, •), we have the following equality: 

Ae(C*) n 1-1 

So let g be so that the above minimum is attained. Without loss of generality, we may scale 
the g l so that Ai = • • • = A n = 1. In that case, 

Wf-gT ^ EUf-gT 

Fill 2 ' 



We are then in the setting of IBCSS98] , pp. 248-250] where we identify our linear / with a 
normally distributed (n+ 1) x n complex matrix. The right-hand side in the above inequality 
is then precisely the left-hand term in ||BCSS98| , Rem. 2, Pg. 250]. Therefore, using the 
notation of ||BCSS98| , Prop. 4], d P (f, S( p , g )) > ^f(/; So it follows that 

z/ Li >, e) = Prob [d P (f, < e] < Prob [d F (f, H x ) < e] 

and the last probability is bounded above by n 3 (n + 1)(#A — 1)(#A — 2)e A via |[BCSS98 
Thm. 6, Pg. 254]. Theorem]]] now follows. ■ 



3.2 The Proof of Theorem @ 

Using |[BCSS98| , Theorem 5 p. 243] (or Proposition ||, Pg. [25] below), we deduce that the 
average number of complex roots is: 



Avg 



n 



-iirn 2 /2 



(2tt) m > 



- det (DG M DGl q) ) 



'(p,q)eU J f&Ffaq) 

By Lemma ^ we can replace the inner integral by a 2n-form valued integral: 



Avg 



_l)n(n-l)/2 



(p,q)EU Jf£F (p>q) 



A 



-n/ii 2 /2 



(27T 



I Mi 



-f ■ ( D VA t ) iP , q) dpA 



Af-(Dv Ai ) iPi _ q) dq . 



Since the image of Dv Ai is precisely (FaJo^) C Ta^ one can add n extra variables 
corresponding to the directions vaXv + qV~1) without changing the integral: we write 
^Ai = ^Ai,(p,q) x ^vaXp + qV-^)- Since (/* + + qV~^)) Dv A . is equal to f l Dv Al , the 

average number of roots is indeed: 



Avg 



\n(n-l)/2 



3 -nrn 2 /2 



/\ (27r)Mi+l fi ' {DVAi)MdpA 



Af -{Dv Ai ) {v> _ q) dq . 
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In the integral above, all the terms that are multiple of f^ffa f° r some a ^ (3 will cancel 
out. Therefore, 



Av g = (-i r cn-i)/2 1 i ^^^j2\mDv A x, q) dP^ 



A (Dv A X,- q) dq 



Now, we apply the integral formula: 



HNP/2 



Fi 



(27T) 



\Xi\ 



-NI 2 /2 
2tt 



to obtain: 



Avg 



(_l)n(n-l)/2 



7T" 



(P,?)e£/ a 



According to formulae |2.1.3| and [2.1.4| , the integrand is just 2 n f\ ui Ai , and thus 
Avg 



(_l)«(«-l)/2 , A „, 



7T' 



/\u Ai = - n I dT« 



3.3 The Proof of Theorem 



Let (p, q) G T n and let / G ^F(p,q)- Without loss of generality, we can assume that / is scaled 
so that for all i, \\f l \\ = 1. 

Let Sf G Ffaq) be such that / + Sf is singular at (p, g), and assume that ll^/*l| 2 * s 
minimal. Then, due to the scaling we choose, 



Pi||2 



Since f + Sf is singular, there is a vector u ^ such that 

'(/ 1 + 5/ 1 )-(^ 1 ) (Pi9) " 



and hence 



{p + 5n-{Dv An ) (m) 

(r+*r)-(i^o<M) 



it = 



it = . 



This means that 



Z 1 ■ Dv Al u = -Sf 1 ■ Dv Al u 



f n -Dv An u = -Sf n -Dv An u 
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Let D(f) denote the matrix 



D{f) 



def 



Given v = D(f) u, we obtain: 



f l ■ (Dv Al ) M 

r ■ (Dv An ) M 

-5f ■ Dv Al D{fY 1 v 



(3.3.1) 



>j n = -5f n ■ Dv An D(f)- 1 v 

We can then scale u and v , such that \\v\\ = 1. 

Claim 1 Under the assumptions above, 5f % is colinear to (Dv Ai D(f)~ 1 v) H . 

Proof: Assume that 5f l = g + h, with g colinear and h orthogonal to (Dv Ai D(f) v) . As 
the image of Dv Ai is orthogonal to v Ai , g is orthogonal to v A ., so ev(g\ (p, q)) — and hence 
ev(h\ (p,q)) = 0. We can therefore replace Sf l by g without compromising equality (|3.3.1| ). 
Since \\Sf\\ was minimal, this implies h = 0. M 

We obtain now an explicit expression for 8f l in terms of v. 



6f 



(Dv Ai D(f)- l v 
\\Dv Ai D{f)-^v\ 



Therefore, 



\\Dv Ai D(f)-iv\ 
So we have proved the following result: 



Lemma 4 Fix v so that \\v\\ = 1 and let Sf E ^(p,q) be such that equation ( \3. 3. 1\ ) holds and 
1 1 5f 1 1 is minimal. Then, 

Lemma ^ provides an immediate lower bound for = v/$Z ll^/ l ll 2: Since 

\\sf\\ > 



max., \\D(f) 1 v\\ Aj 



we can use \\v\\ = 1 to deduce that 



i 2 > 



> 



maxj \\D(f) 1 v\\ A] m&Xj\\D(f) 1 || A 
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Also, for any v with \\v\\ = 1, we can choose 5f minimal so that equation ( |3.3.1| ) applies. 
Using Lemma |J we obtain: 

ll*/*ll< — 



Hence 



'En*/' 



i 2 < 



minj WDif^vllAj 
Since this is true for any v, and \\Sf\\ is minimal for all v, we have 

1 



I 2 < 



max||„|| = i minj 1 \\ Aj 
and this proves Theorem [|. 

3.4 The Idea Behind the Proof of Theorem [5] 

The proof of Theorem [5] is long. We first sketch the idea of the proof. Recall that Fop q) is 
the set of all / G T such that ev(f;p + = 0, and that is the restriction of the 

discriminant to the fiber Fop^)'- 

S( P g) = f {/ G ^F{p, q ) '■ does not have full rank} . 



The space T is endowed with a Gaussian probability measure, with volume element 

e -||/H 2 /2 

where dT is the usual volume form in T = {Ta x -, (•, x • • • x (^4 n) (•, ■},*„) an d ||/|| 2 = 
^ H/ 1 !!^.. For £7 a set in T n , we defined earlier (in the statement of Theorem |5|) the quantity: 

u A (U,e) = Prob[Ai(/,E0 > ^T 1 ] = Prob[3(p,g) G £7 : d ¥ (f,Z M ) < e) . 

The naive idea for bounding i/ A (Z7, e) is as follows: Let V(e) = {(/, (p, ? )) G F x U : 
ev if](PyQ)) — and dp(/, E(p )9 )) < e}. We also define 7r : V(e) — > T as the canonical 

projection mapping J 7 x U to J 7 , and set #v(e)(f) = #{(p,q) e J7 : (f,(p,q)) G V(e)}. 
Then, 



X7r(V(e))(/) 



"72 



(2vr)E 



< 



172 



(2tt)E 
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dT 



with equality in the linear case. 



Now we apply the coarea formula IBCSS98] , Theorem 5 p. 243] to obtain: 



ap(/,S( p>g ))<e 



(p,q)£UcT n 



where dVr« stands for Lebesgue measure in T n . Again, in the linear case, we have 
equality. 

We already know from Lemma |3] that 

n 

l/NJ{f; (p, q)) = /\f- (Dv Ai ) M d P A f ■ (Dv Ai ) M dq . 
i=i 

We should focus now on the inner integral. In each coordinate space T Av we can introduce 
a new ortho normal system of coordinates (depending on (p,q)) by decomposing: 

/' /; • ./;; • /; , 

where /* is the component colinear to v A ., fl is the projection of f l to (range Dv Ai ) H , 
and is orthogonal to f* and 

Of course, f % G (^ r A i )(p,, 3 ) if and only if /' = 0. 
Also, 

n 

A /* • ( Dv A t )( P , q )dp A f ■ (Dv Ai ) M dq = 

n 

= A ' ( Dv Ai)(p,g) d P A fl ■ ( D VA t )(p,q)dq . 



i=l 



i=l 



It is an elementary fact that 
It follows that for / e F{p, q )'- 
with equality in the linear case. Hence, we obtain: 



u A (U,e)< I I fer(pq) l/\f:-(Dv At ) M dpAf:-(Dv At ) M dq 

R -n/i+/;ii 2 /2 



with equality in the linear case. We can integrate the ^(Mj — n — 1) variables / n to obtain: 
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Proposition 3 

v A (U,e)< [ [ (j\fl-(Dv At ) M dpAfl.(Dv A J M dq 



'( P ,<?)eC/cT" _ v 

ap(jn,^( Pl9 ))<£ 



v«=l 

e" 



ll/^l| 2 /2 



n(n+l) 



(2 7 r) n ( n + 1 ) 
uni/i equality in the linear case. ■ 

3.5 The Proof of Theorem |5] 

The domain of integration in Proposition ^ makes integration extremely difficult. In order 
to estimate the inner integral, we will need to perform a change of coordinates. 

Unfortunately, the Gaussian in Proposition ^ makes that change of coordinates extremely 
hard, and we will have to restate Proposition § in terms of integrals over a product of 
projective spaces. 

The domain of integration will be P n_1 x • • • x P n_1 . Translating an integral in terms of 
Gaussians to an integral in terms of projective spaces is not immediate, and we will use the 
following elementary fact about Gaussians: 

Lemma 5 Let (p : C n —>■ R be C* -invariant (in the sense of the usual scaling action). Then 
we can also interpret ip as a function from P n_1 into R ; and: 

1 f f e -M\ 2 n 

<p(x)d[x] = / (p(x) . dx , 



Vol(P™+ 1 )i N6P , l - 1 " v ' 1 1 A ec „" v ' (2tt)« 

where, respectively, the natural volume forms on P™ -1 and C™ are understood for each inte- 
gral. 

Now the integrand in Proposition ^| is not C*-invariant. This is why we will need the 
following formula: 

Lemma 6 Under the hypotheses of Lemma 

If If e -IMI 2 /2 

<p(x)d[x] = — / ||a;|| 2 93(a;)— — - — dx . 



Vol(P™+ 1 )i Wepil - 1 ^ v L 1 2nA et> " rw (2tt 

where, respectively, the natural volume forms on P ra_1 and C n are understood for each inte- 
gral. 



Proof: 

/ ||x|rV(x) ". dx = I I [r[ 2n+ V(Q) " /n x drde 



e -IMI 2 /2 r r«> e"l r l 2 / 2 

12,^^ Arr. — I / U|2n+1 - 



xec n 



(2tt)» Jecs^Jr=o (2tt)» 
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r ( e -M72 °° roo f -\r\ 2 /2 \ 





/• e -|NI| 2 /2 

= 2n (f(x)—, — r^-dx 
Jxec n (27r) n 

■ 

We can now introduce the notation: 

n 

WEDGE A (/ n ) d ^ /\ ——Jl ■ (Dv Ai ) M dp A £ • {Dv Ai ) (p , 
i=i ii/eI1 

This function is invariant under the (C*) n -action \~k fx '■ /[ 1 — > (Ai/j, • • • , A n /™). 
We adopt the following conventions: JF E c is the space spanned by coordinates / n and 
POFh) is its quotient by (C*) n . 

We apply n times Lemma || and obtain: 

Proposition 4 Let VOL = f Vol(P n_1 ) n . T/jen, 

< / / f ePf ^ , WEDGE A (/ E ) dF(F E ) dF T n 



rfp(/l>S( p>9 ))<£ 



and in the linear case, 



Now we introduce the following change of coordinates. Let L G GL(n) be such that the 
minimum in Definition [5| p. [| is attained: 

tp : P 71 " 1 x • • ■ x P™- 1 -> P"- 1 x ■ • • x P™- 1 

/n i-> fli = f <p(/i) 3 such 

that g-j = f* ■ Dv Al L . 

Without loss of generality, we scale L such that det L = 1. The following property follows 
from the definition of WEDGE: 

WEDGE A (/ n ) = WEDGE Lin (^ E ) TT (3.5.1) 

Assume now that dp(/ n , < e. Then there is 5f G JF E , such that / + 5f G S^ 11 ^ and 

1 1^/11 < e (assuming the scaling ||/*|| = 1 for all i). 



Setting g % = ip(f z ) and Sg = f(g), we obtain that g + 5g G S L 



0,9)' 



^,sg n 9) )< 



IN 
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At each value of i, 

WW / \\sf\ 



< 



\9i\\ - 



where k denotes Wilkinson's condition number of the linear operator Dj^p 1 . This is precisely 
n(DvAiL). Thus, 

d P (g, S^ n ?) ) < emaxn(Dv At L) = max ^/n(uj Ai ) 



Thus, an ^-neighborhood of ^f pq ^ is mapped into a y/ k\jE neighborhood of S^ 1 ^. 
We use this property and equation ( ^.5.1| ) to bound: 



(2n) n 



- A (U,e) < ^L- I L^!,...^-! WEDGE- 

A 112 



VOL ( va ^ucT- J n&n ~, x '" xl 



/I/ 



■ n w 1 ^^ 1 ' 2 rf(pn_l x • • • x pn_l) dVrn (3 - 5 - 2) 

where Jg^f' 1 is the Jacobian of <^ _1 at g E . 
Remark 3 Considering each Dv Al as a map from C n into C n , the Jacobian is: 



\f Hg* Y\\ n ^.h n „. ^1/2 

Ils5ll : 

H^e u>z// not nse this value in the sequel, o 



'n*- 1 = W V (4etDv«Dv A ,y 

1=1 



In order to simplify the expressions for the bound on v (U, e), it is convenient to introduce 
the following notation: 

def (2n) " F Lin T n ^ (Pn " X X • ' ' X ^ 
dP " VOL WEDGE M „! K in )An 

n II * 1 1 2 

* - ys |j "' 

def 

X5 " X {9:*(3,E^ n 9) )<5} 

Now equation ( p.5.2| ) becomes: 



zr 4 ^) < n! / (w Lin ) A " / x^ £ (a n ) (3.5.3) 

Lemma 7 Let (p, g) be fixed. Then P n_1 x • • ■ x p n_1 together with density function dP, is 
a probability space. 
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Proof: The expected number of roots in U for a linear system is 



An 
Lin 

(p,q)eU Jg I eP n - 1 x---xP n - 1 



n\ I / dP 



which is also n\ jjjU^- This holds for all U, hence the volume forms are the same and 

dP = 1M 



3neP n - 1 x---xP n - 1 



This allows us to interpret the inner integral of equation ( |3.5.3| ) as the expected value of 



a product. This is less than the product of the expected values, and: 
v A (U,e) <n\ f K in )A« (f dP H( 9ll )) • 

■ ( dp x^ue(g*) J 

\Jg I eP' l - 1 x---xP"- 1 / 



Because generic (square) systems of linear equations have exactly one root, we can also 

(to 1 



consider U as a probability space, with probability measure J n , . ralo;^" . Therefore, we 



can bound: 



1 



» A (U,e) < ^0777777 / n\{u hm )^ / dP H^] 

VOl {(J) \J(p,q)£U J g E eP"- 1 x-xP"- 1 



n\(u U n) An / dP X ^Ue{g.) 
{p,q)£U J 3neP n - 1 x---xP n - 1 

The first parenthetical expression is Vo\ A (U), the volume of U with respect to the 
toric volume form associated to A = {A\, • • • ,A n ). The second parenthetical expression 
is ^ Lin (y / K7j r e, U). This concludes the proof of Theorem || ■ 



3.6 The Proof of Theorem |3] 

As in the complex case (Theorem |2|), the expected number of roots can be computed by 
applying the co-area formula: 



AVG= / / TT- j^detiDG DG 11 )- 1 / 2 

JpeuJfe^fJl A/27T ' 



" -n/ii 2 /2 

* JL. I 1 

IpeuJferfi-x V27T 

Now there are three big diferences. The set U is in M n instead of T ra , the space J-^ 
contains only real polynomials (and therefore has half the dimension), and we are integrating 
the square root of 1/ det (DG DG H ). 

Since we do not know in general how to integrate such a square root, we bound the inner 
integral as follows. We consider the real Hilbert space of functions integrable in J 7 ^ endowed 
with Gaussian probability measure. The inner product in this space is: 
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where dV is Lebesgue volume. If 1 denotes the constant function equal to 1, we interpret 
AVG — [ (2n)- n / 2 (det(DG DG H )- 1/2 ,1) . 

JpEU 



i P eu 

Hence Cauchy-Schwartz inequality implies: 



AVG < [ (2n)-^ 2 \\ det(DG DG H y 1/2 



'pec 

By construction, ||1|| = 1, and we are left with: 



r r p-ll/^ll 2 /2 

AVG< / (2vr)-"/ 2 / JJ_^^det(DGDGH)-i . 



As in the complex case, we add extra n variables: 



r r 71 ~ ll/ i ll 2 /2 
AVG<(2n)- n / 2 / Y\ ~~=m~ det(DG DG H ) _1 , 



and we interpret det(_DG DG H ) 1 in terms of a wedge. Since 

2 e-iNi 2 / 2 r 2 e~y 2 / 2 r e~y 2 /' 

^/2~K M JyeR J y eR V2tt 



we obtain: 

AVG < (27r)- n / 2 f VridT^ = (2n)- n/2 [ VridT^ . 

Jp&U Jp& 

Now we would like to use Cauchy-Schwartz again. This time, the inner product is defined 

as: 

{(p,rf>) = / ip(p)ip(p)dV . 
J P eu 

Hence, 

AVG < (2n)- n/2 (n\dT n ,l) < (2n)- n/2 \\n\dT n \\ ||1|| . 
This time, ||1|| 2 = X(U), so we bound: 



AVG < (2it'y n/2 ^MU)^ J n\dT n 

< (An 2 )- n/2 ^/MU)J [ n\dT« 

Y J(p,q)eT",p£U 
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3.7 The Proof of Theorem © 

Let e > 0. As in the mixed case, we define: 



def 



Prob /e ^ [ t i(f;U)>s- 1 ] 

Prob /e ^r [Bp e U : ev(f;p) = and d P (f, £ p ) < e] 



where now f/el™ 

Let V(e) = {(f,p) £ F K x U : ev(f;p) = and dp(/, £ p ) < e}. We also define 
7r : V(e) — > P(JF) to be the canonical projection mapping Fr x U to Fr and set #vYe)(/) = 
#{pet/:(/,p)ey(e)}. Then, 



< 



< 



^_ FM . X»r(V( g ))(/) ^ 

e-Eill^HVa 

/• e -EJI/il 2 /2 ! r 



As before, we change coordinates in each fiber of J- A by 

f = fi + f* + fm 

with /* colinear to 1/4, (/j) T in the range of Dva, and fl othogonal to /* and f*. This 
coordinate system is dependent on p + T. 

In the new coordinate system, formula |2.3.1| splits as follows: 

det (DG ip) DGf p) y 1/2 dVrn = 



det 



det 



(/l )n 

(/n )n 
(/i )n 

(/i )n 



det 



det Dvf -Dv A 



(Dv A 



n 



n 



The integral E(U) of a/ det DvaDv a is the expected number of real roots on U, therefore 

e -£dl/l+/^l72 



i (U,s)<E(U) 



rfp(/l+/M,S p )<£ 



det 



• • • (A)n 
(/iT)l • • • 



23 



In the new system of coordinates, S p is denned by the equation: 



det 



(/n)l ••• (/i)n 
{f£)l ■ ■ ■ (fx)n 



Since HA + /.H > ||M 

This implies: 

>u(U,e)<E(U) 



dp(U + U, S p ) < e =>• d P (f JL , Ep) < e 



e -EJl4+/ill 2 /2 



d P (/n,[det=0])<e VZ/I 



det 



"a, 1 : 



(/■ n 



(/■" 



We can integrate the Mj — n — 1) variables / M to obtain: 

e-Eill/lll 2 ^ 



l (U,e) = E(U). 

d P (/n,[det=0])<e V -" 



det / n | 



This is E(U) times the probability u(n,e) for the linear case. 
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A The Coarea Formula 



Here we give a short proof of the coarea formula, in a version suitable to the setting of this 
paper. This means we take all manifolds and functions smooth and avoid measure theory as 
much as possible. 

Proposition 5 

1. Let X be a smooth Riemann manifold, of dimension M and volume form \dX\. 

2. Let Y be a smooth Riemann manifold, of dimension n and volume form \dY\. 

3. Let U be an open set of X, and F : U — > Y be a smooth map, such that DF X is 
surjective for all x in U . 

4- Let ip : X — > M + be a smooth function with compact support contained in U. 

Then for almost all z G F[U), V z == F~ 1 (z) is a smooth Riemann manifold, and 

cp{x)NJ(F;x)\dX\ = I [ cp(x)\dV z \\dY\ 

JzeY Jx&V z 

where \dV z \ is the volume element of V z and NJ(F,x) = a/ det DF^DF X is the product of 
the singular values of DF X . ■ 

By the implicit function theorem, whenever V z is non-empty, it is a smooth (N — n)- 
dimensional Riemann submanifold of X. By the same reason, V := {(z, x) : x G V z } is also 
a smooth manifold. 

Let f] be the following iV-form restricted to V: 

r] = dY A dV z . 

This is not the volume form of V. The proof of Proposition [5] is divided into two steps: 
Lemma 8 

j ip(x)\r]\ = [ <p(x)NJ(F;x)\dX\. 
Jv Jx 

Lemma 9 

/ V (x)\ V \= [ [ V (x)\dV z \\dY\ . 

Jv JzeY JxeV z 

Proof of Lemma §|: We parametrize: 



Then, 



i>\ x -> v 

x t— * (F(x),x) 



<p(x)\r}\ = / (ipoip)(x)\ip*T]\ . 
\ Jx 
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We can choose an orthonormal basis tti, • • • ,% of T X X such that u n +i, ■ • • ,um € 
kerDi 7 ^. Then, 

(DF x u h Ui) i = 1, • • • ,ra 
(0,«i) i = n+l,---,M ' 



Thus, 



|^*r/(wi, • • • ,mm)| = ^(D^wi, • • • ,DV^m)| 

= \dY{DF x u x , ■ ■ ■ , DF x u n ) \ \dV z (u n+1 , ■ ■ • ,u M ) 

= | det DF x \ kevDF ±\ 

= NJ(F,x) 



and hence 



ip(x)\r]\ = / <p(x)NJ(F;x)\dX\. 



Proof of Lemma |9|: We will prove this Lemma locally, and this implies the full Lemma 
through a standard argument (partitions of unity in a compact neighborhood of the support 
of (p). 

Let x ,z be fixed. A small enough neighborhood of (x ,z ) C V Zo admits a fibration 
over V ZQ by planes orthogonal to kerDF Xo . 
We parametrize: 

9: YxV Z0 ^ V 

(z,x) ^ (z,p(x,z)) ' 

where p(x,z) is the solution of F(p) = z in the fiber passing through (zo,x). Remark 
that 6*dY = dY, and 6*dV z = p*DV z . Therefore, 

6*(dY A dV z ) = dY A (p*dV z ) . 

Also, if one fixes z, then p is a parametrization V ZQ — > V z . We have: 

ip(x)\r]\ = / ip(p(x,z))\e*rj\ 



v 



<p( P (x,z)\p*dv z \ ] k/r| 



The proposition below is essentially Theorem 3 p. 240 of |BCSS98|. However, we do not 



require our manifolds to be compact. We assume all maps and manifolds are smooth, so 
that we can apply proposition [5]. 

Proposition 6 

1. Let X be a smooth M- dimensional manifold with volume element \dX\. 
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2. Let Y be a smooth n- dimensional manifold with volume element \dY\. 

3. Let V be a smooth M -dimensional submanifold of X x Y , and let 7Ti : V — > X and 
7r 2 : V — > Y be the canonical projections from X x Y to its factors. 

4- Let E' be the set of critical points of tci, we assume that £' has measure zero and that 
£' zs a manifold. 

5. We assume that tt 2 is regular (all points in tt 2 (V) are regular values). 

6. For any open set U C V , for any x G X, we write: #u(%) = #{ 7r r 1 ( a; ) H f/}. We 
assume that f xeX #v(x)\dX\ is finite. 

Then, for any open set U C V , 

[ #u(x)\dX\= I [ ^==L==\dV z \\dY\ 
Jx^iu) Jz& J^Ijj \/det DG X DG* 

where G is the implicit function for (x, G(x)) & V in a neighborhood of (x, z) G V \ E'. ■ 

Proof: Every {x, z) G U \ £' admits an open neighborhood such that 7Ti restricted to that 
neighborhood is a diffeomorphism. This defines an open covering of U \ E'. Since £7 \ E' is 
locally compact, we can take a countable subcovering and define a partition of unity (<^a) AeA 
subordinated to that subcovering. 

Also, if we fix a value of z, then (<^a)aga becomes a partition of unity for 7ri(7Tf (V z ) PI £/). 
Therefore, 



/ 



#C7(x)|dX| = X) / ^A(x,^)|dX| 



V / / ^fo*) u x \ 

^ I - ' ;)ZeSup p^ NJ(G,x) 



AeA 



\- / vA(-r-c) |/v| 



zey Xs= a ^i',2eSuppi/3 A NJ(G,x 



AeA " a,*teoupp^ 

1 



NJ(G,x] 



\dX\ 



where the second equality uses Proposition |5| with tp = (f\/NJ. Since NJ= ^/det DG X DG*, 
we are done. ■ 
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